

%% parameters
load SMM
fprintf('verif poids=1');
p0 =[ p01_SMM;p02_SMM;p03_SMM;p04_SMM];
mu_c = [mu_c1_SMM; mu_c2_SMM; mu_c3_SMM; mu_c4_SMM];
%sig_c = sig_c_SMM;
sig_eps_a = [sig_eps_a1_SMM; sig_eps_a2_SMM; sig_eps_a3_SMM; sig_eps_a4_SMM];
%rho_a = [rho_a1_SMM; rho_a2_SMM; rho_a3_SMM; rho_a4_SMM];


    [f1, avfracup1,    med1,    interq1, kur1, f2, avfracup2,    med2,    interq2, kur2,f3, avfracup3,    med3,    interq3, kur3,f4, avfracup4,    med4,    interq4, kur4]=geNCalvoPlus_SMM4(p0, mu_c, sig_eps_a, weight_j);
    
    
    
    